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Abstract. 

We present the preliminary results of the analysis to search for inspiraling 
compact binaries using TAMA300 DT8 data which was taken during 2003. We 
compare the quality and the stability of the data with that taken during DT6 in 
2001. We find that the DT8 data has better quality and stability than the DT6 
data. 



PACS numbers: 95.85.Sz, 04.80.Nn, 07.05Kf, 95.55Ym 

TAMA300 fX] is an interferometric gravitational wave detector with 300 m baseline 
length located at Mitaka campus of the National Astronomical Observatory of Japan 
in Tokyo (35.68° N, 139.54° 'E). TAMA300 began to operate in 1999. It had performed 
observations for 8 times by July 2003. In particular, during the period from 1 August 
to 20 September 2001, TAMA300 performed an observation, which is called Data 
Taking 6 (DT6), and about 1039 hours of data were taken. The best sensitivity 
was about 5 x 10~ 21 /\/Hz around 800Hz. More recently a longer observation was 
performed by TAMA300 during 14 February to 15 April. 2003, and about 1163 
hours of data were taken. This observation is called Data Taking 8 (DT8). The 
best sensitivity was about 3 x 10~ 21 /yTIz around 1kHz. 

We consider gravitational waves from inspiraling compact binaries, consisting of 
neutron stars or black holes. Since their wave forms are well known by post-Newtonian 
approximation of general relativity we can use the matched filtering. In this paper, 
we report some results applied to DT6 and DT8 data, and compare the stability of 
both observations. The data used for the analysis are that taken during 1 August to 
20 September, 2001 (DT6), and during 14 February to 15 April, 2003 (DT8). The 
length of data is 1039 hours and 1163 hours for DT6 and DT8 respectively. These are 
the data after removing any unlocked period, the period when some adjustments were 
made, and the locked period continued only less than 10 minutes. 

We assume that the time sequential data of the detector output s(t) consists of 
a signal plus noise n(t). We denote the one-sided power spectrum density of noise by 

Sn(f). 
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The gravitational wave strain amplitude h(t) is calculated by combining two 
independent modes of the gravitational wave and antenna pattern of interferometer 
as 

h(t) = A[h c (t - t c ) cos 0o + h s (t - t c ) sin0 o ], (1) 

where h c (t) and h s (t) are two independent templates whose phase are different to 
7r/2 from each other and A is amplitude. The parameters which characterize the 
templates are the coalescence time t c , the chirp mass A4(= Mrf/ b ) (M — m\ + m?), 
non-dimensional reduced mass r]{= mim2/M 2 ) and the phase <fio. The details of the 
wave form can be found in [2J. We use the restricted post-Newtonian wave forms in 
which the phase evolution is calculated to 2.5 post-Newtonian order, but the amplitude 
evolution contains only the lowest Newtonian quadrupole contribution. The effects of 
spin angular momentum are not taken into account here. 

The two wave forms, h c and h s , are transformed into the Fourier domain by the 
stationary phase approximation. We use templates for matched filtering defined in the 
Fourier domain as h(f) = h c (f) cos^o + h s (f) sin0O: where h(f) denotes the Fourier 
transform of h(t). Here, h c and h s are normalized as (h c \h c ) = (h s \h s ) = 1, where 

<*-/>™ 

and * means the operation to take the complex conjugate. The filtered output is then 
defined by 

p(t Ci M,n,(h) = (s\h). (3) 

In the matched filtering, we look for the maximum of p over the parameter t c , M, rj, 
and 0o- In Eq. J3J), we can take the maximization over </>o analytically which gives 



p{t c ,m 1 ,m 2 ) = y/(s\h c ) 2 + {s\h s ) 2 , (4) 

We see that p 2 has an expectation value 2 in the presence of only Gaussian noise. 
Thus, we can define the normalized signal-to- noise ratio by SNR = p/^/2. 

In the TAMA300 analysis, we have found that the noise contained a large amount 
of non-stationary and non-Gaussian noise 3 . In order to remove the influence of such 
noise, we introduce a x 2 test [5]. First we divide each template into n mutually 
independent pieces in the frequency domain, chosen so that the expected contribution 
to p from each frequency band is approximately equal, as 

Ks(f) = h^{f) + U 2 l(f) + ■■■ + %)(/). (5) 

The x 2 is defined by 



X 2 = 

i=l 

where 



This quantity must satisfy the x 2 -statistics with 2n — 2 degrees of freedom, as long 
as the data consists of Gaussian noise plus inspiraling signals. For convenience, we 
use the reduced x _sc i uar e defined by x 2 /(2n — 2). In this paper, we chose n = 16 
(hereafter, x 2 means the reduced x -sc iuare). 

The value of x 2 1S independent of the amplitude of the signal as long as the 
template and the signal have an identical wave form. However, in reality, since the 
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template and the signal have different value of parameters because of the discrete time 
step and discrete mass parameters we search, the value of x 2 becomes larger when 
the amplitude of signal becomes larger. In such situation, if we reject events simply 
by the value of x 2 , we may lose real events with large amplitude. In order to avoid 
this, we have been introducing a statistic, pj \fy? , in the TAMA300 data analysis JBj. 
The event selection using pi \[y? is more efficient to detect events with large p (say, 
p > 20) than the usual x 2 selection. Further, by checking the Galactic event detection 
efficiency, we have found that the detection efficiency increases about 20 % by the 
pj \[y? selection compared to the x 2 selection in the case if the threshold is x 2 < 1.5. 

The mass parameters considered in this paper is 1.0M Q < mi,m,2 < 2.0M Q , 
which is a typical mass region of neutron stars. In this region, a discrete mass 
parameter space is determined so that the maximum loss of SNR due to discretization 
becomes less than 3%. The parametrization of mass parameters are done based on 



The mass parameter space depends on the power spectrum of noise. Averaged 
power spectrum of noise for each continuously locked segment was used to define 
the mass parameter space in each segments. This parameter space is not changed 
within the continuously locked segment. However, in order to take into account of the 
variation of the noise power spectrum with time, we use a different mass parameter 
space for different locked segments. The mass parameter space turned out to contain 
about 200 ~ 1000 templates for the DT6 data, and about 200 ~ 800 templates for the 
DT8 data (Fig.l). The typical value of the number of template is about 700 for DT6, 
and 600 for DT8. The variation of the number of template is not due to the variation 
of absolute amplitude of the power spectrum, but due to the variation of the shape of 
the power spectrum. Thus we can see an aspect of variation of the shape of the power 
spectrum from Fig. 1 . It is found that DT8 data are more stable than DT6 data with 
respect to the number of templates, which probably means the stability of DT8 data 
with respect to the shape of the power spectrum. 

Next we examine the variation of total noise power. Here, we define 



power is suited to examine the variation of the noise power spectrum used in the 
matched filtering. The value of p is evaluated for each data with length 1.1 minutes. 
Here, we only consider the locked data. In Fig. 2, we plot the number of occurrence 
as a function of the value p. The mean value and standard deviation are 3.2 x 10 -19 , 
1.5 x 10" 19 for DT6 and 1.9 x 10" 19 , 7.5 x 10" 20 for DT8 respectively. It is found 
that the stability of DT8 data is better than that of DT6 data with respect to the 
amplitude of the noise power spectrum. 

Next we present results of matched filtering analysis including preliminary results 
of DT8 analysis. After calculating Eq.Q for each mass parameter in the mass 
parameter space, we search for t c and masses which realize the maximum value of 
p in each interval of the coalescence time with length At c — 25.6 msec. For those t c 
and masses, the value of the reduced x-square are calculated. In Fig. 3, we plot the 
number of events as a function of p and yx 2 . In these figure, only events with p > 7 
are plotted. We find from the distribution of p that the DT8 result has longer tail 
than DT6 data. Such events with large p must be due to the non-Gaussian noise since 
the value of yx 2 °f them are also very large. The distribution of yx 2 is a l so more 




-1/2 



(8) 



Since this value is related to the normalization constant of the templates, this noise 
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spread in DT8 case than in DT6 case. Thus, in terms of the distribution of p and 
y/X^, the stability of DT8 does not seem to be better than DT6. However, most of 
events with large p in DT8 have larger yx than that of DT6. Thus, if we plot the 
number of events as a function of pj y/x?, we have a shorter tail in DT8 than in DT6. 
This can be seen in Fig. 4. If we determine threshold of pj yf)d f° r a given false alarm 
rate, the threshold of DT8 becomes lower than that of DT6. This fact will help us 
to improve the detection efficiency of the Galactic events in DT8 further. These facts 
also suggest that the DT8 data have much better quality than the DT6 data in the 
sense that it is easier to discriminate the non-Gaussian noise events. 

After the observation of DT6, the configuration of TAMA300 has been changed 
very much. One of the major changes is the installation of the power recycling system. 
Thus, it will be important and interesting to investigate the origin of the different 
property of data between DT6 and DT8. We will work on such analysis in the future. 

Much more details of the analysis using the data of TAMA300 DT6 and DT8 will 
be discussed elsewhere UJ. 
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Figure 1. The frequency distribution of the number of templates used in DT6 
and DT8 analysis are plotted 
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Figure 2. Noise power p for each locked data with 1.1 minutes. The left and 
right panel shows the average of noise power for DT6 and DT8 respectively. 
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Figure 3. The distribution of the number of events are plotted as a function of 
P and v 7 ^". 
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Figure 4. The number of events as a function of p/\/x 2 - The left and right 
panel shows histograms for DT6 and DT8 respectively. 



